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' _ ^ Abstract. We describe an algorithm that recognizes some (perhaps all) in- 

trinsically knotted (IK) graphs, and can help find knotless embeddings for 
graphs that are not IK. The algorithm, implemented as a Mathematica pro- 
gram, has already been used in ^ to greatly expand the list of known minor 
minimal IK graphs, and to find knotless embeddings for some graphs that had 
previously resisted attempts to classify them as IK or non-IK. 

in 

H 
O 

All graphs in the paper arc finite. A graph is intrinsically knotted (IK) if 
every embedding of it in contains a nontrivial knot. In the early 1980s Conway 
and Gordon [T] showed that Kj, the complete graph on seven vertices, is IK. A 
^ graph _ff is a minor of another graph G if H can be obtained from a subgraph 

I— I of G by contracting some edges. A graph G with a given property is said to be 

^ minor minimal with respect to that property if no proper minor of G has that 

^ property. Robertson and Seymour's Graph Minor Theorem [9] says that in any 

infinite set of graphs, at least one is a minor of another. It follows that for any 
property whatsoever, there are only finitely many graphs that are minor minimal 
with respect to that property. In particular, there are only finitely many minor 
' I minimal IK graphs. It is well known and easy to show that G is IK iff G contains 

0^ a minor minimal IK graph as a minor. Hence, if we had a list of all minor minimal 

IK graphs, then we could determine whether or not any given graph G is IK by 
^ I testing whether G contains one of the graphs in this list as a minor. But, for about 

• . the past three decades, finding such a list has remained an open problem. Our 

algorithm, which does not rely on having such a list available, can recognize some 
(perhaps all) IK graphs. 

Prior to this algorithm, finding new minor minimal IK graphs was quite difficult. 
^ Initially, the only known minor minimal IK graphs were ify and the thirteen graphs 

obtained from it by VY moves ([T], [S]). About two decades later the list expanded 
to include ^^3,3,1,1 plus the 25 graphs obtained from it by VY moves ([5], [S]). And 
a few years later Foisy [S] discovered one more minor minimal IK graph. These 
41 graphs were the only minor minimal IK graphs known before this algorithm. 
With the aid of our algorithm, Goldberg, Mattman, and Naimi [6] have found over 
1600 new IK graphs, over 200 of which have been proven to be minor minimal IK. 
The algorithm has also helped us find knotless embeddings for some graphs that 
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had previously resisted attempts to determine whether or not they are IK. For the 
graphs we worked on, our program typicahy took from less than a minute to less 
than a day to run on an average personal computer. A large portion of the time is 
usually spent on just finding all the cycles in the given graph; and a good predictor 
of computation time seems to be the number of edges in the graph. We found that 
for some graphs with about 30 or more edges, the program could take days (or 
much longer) to run. 

Before describing our algorithm, we need to introduce some more notation and 
terminology. The linking number of two disjoint cycles Ci and C2 in a spatial 
graph (i.e., a graph embedded in S^) is denoted by lk(Ci, C2). If K is a knot in S'^, 
a2{K) denotes the second coefficient of the Conway polynomial of K. The graph 
depicted in Figure [T] is denoted by D4. It has four vertices and eight edges. The 



cycles Ci , • • • , C4 are always assumed to be labeled in the order shown, so that 
Ci n C3 = and C2 n C4 = 0. 



If D4 is embedded in in such a way that lk(Ci, C3) 7^ and lk(C2, C4) ^ 0, 
then we say it is a double-linked D4. If these linking numbers are both nonzero 
in Z2, we say the embedded D4 is double-linked in Z2. The following lemma was 
proved by Taniyama and Yasuhara in Z and by Foisy 2' in Z2. 

Lemma. (Z?4-Lemma, [2], Every double-linked D4 (in "L or 'Li) contains a 

knot K such that 02 (-fC) ^0 (inX or Li, respectively). 

A £'4-less embedding (in Z or Z2) of a graph is an embedding that contains no 
double-linked D4 (in Z or Z2) as a minor (note: a minor H of an embedded graph 
G can in a natural way be given an embedding that is "derived" from the given 
embedding of G). Since any knot K with ai{K) ^ is nontrivial, an immediate 
corollary of the above lemma is: 

Corollary. // a graph has no D4-less embedding (in L or Li), then it is IK. 

Given any graph G, our algorithm searches for a Z?4-less embedding of G in Z 
or Z2. If it finds one, it reports this embedding, which, as we'll explain later, may 
help us find a knotless embedding for G. If the algorithm determines that there is 
no such embedding, we conclude by the above corollary that G is IK. 

Using this algorithm, Goldberg, Mattman, and Naimi 'W have found over 1600 
new IK graphs, most of which seem to be minor minimal IK, and 222 of which were 
"manually" verified in [5] to be minor minimal IK. 




Figure 1. The graph D4. 
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Foisy [4] has asked the following question (which is the converse of the above 
corollary) : 

Question [4]: Does every graph that has a £>4-less embedding (in 
Z or Z2) have a knotless embedding? 

If the answer to this question is affirmative, then in fact our algorithm (without 
the third shortcut described in Section [s]) always decides whether or not a graph 
is IK. In "SI, in every instance that the algorithm found a _D4-less embedding for a 
graph G, whenever the authors tried to manually find a knotless embedding for G, 
they succeeded. 

Let's define a spatial graph F to be an a2-less embedding (in Z or Z2) if for 
every knot K cT, a2(K) = (in Z or Z2). It is then natural to ask: 

(1) Does every graph that has a Z?4-less embedding (in Z or Z2) have an a2-less 
embedding (in Z or Z2)? 

(2) Does every graph that has an a2-less embedding (in Z or Z2) have a knotless 
embedding? 

The algorithm sets up systems of linear equations involving linking numbers of 
certain pairs of cycles in G in terms of variables representing the number of crossings 
between pairs of edges. These equations are such that there is a solution to at least 
one of the systems of equations if and only if G has a £'4-less embedding. The 
number of equations the algorithm needs to solve grows exponentially with the size 
of G. So it was necessary to find various shortcuts to speed up the algorithm. We 
describe the algorithm in Section [2j the shortcuts in Section |3] and the proofs in 
Section [H 

Using similar techniques, we have also written a Mathematica program that 
determines whether or not any given graph is intrinsically linked (IL), i.e., every 
embedding of it in contains a nontrivial link. The algorithm relies on the work 
of Robertson, Seymour, and Thomas lOj, which implies that if a spatial graph 
contains a nontrivial link, then it contains a 2-component link with odd linking 
number. Note that, by [TO. . an algorithm for recognizing IL graphs was already 
known: simply check whether the given graph contains as a minor one of the seven 
graphs in the Petersen Family — which, by [S], can be done in polynomial time 
(though the polynomial has very large coefficients). And in [12 , van der Hoist gives 
polynomial-time algorithms to find an embedding of a graph with all linking num- 
bers zero, given that such an embedding exists. Our algorithm is not polynomial 
time (since the number of cycles can grow exponentially with the number of edges); 
but, in practice, we found it satisfactorily fast and practical for the graphs that we 
tried out: for graphs with 28 or fewer edges, the program usually took from about 
a few seconds to a couple of hours to run on a typical personal computer in 2009. 
More specifically, the computer program was run on at least the following graphs 
listed in [B] and [7]: all graphs in the Kt, K^ ^^ i i, and Eg + e families (188, all 
together); the 25 parentless ("top level") graphs and a few dozen other graphs in 
the Gg^28 family; and the graph Gi4^25 plus a few more in its family. 

2. The Algorithm 

A quad in a graph G is an ordered 4-tuple (Ci, C2, G3, C4) of distinct cycles in 
G such that Ci n G3 = C2 n C4 = 0, and for each i G {1, 2, 3, 4}, either Q n C^+i 
is nonempty and connected (where an index of 4 + 1 is to be interpreted as 1), or 



4 



JONATHAN MILLER AND RAMIN NAIMI 



Ci n Ci+i = and there exists a path pi in G with one endpoint in Ci and the other 
in Ci-|_i, satisfying the following conditions: 

(1) Each path pi is disjoint from Cj for j ^ {i, z + 1}. 

(2) The interior of each pi is disjoint from Ci and C^+i. 

(3) Any two distinct paths pi and are disjoint. 

Furthermore, for each i where Ci H C^+i is nonempty and connected, we let pi 
consist of an arbitrary vertex in Cif]Ci+i. (We do this so that pi is defined for every 
i, allowing for more uniformity in notation). We say Pi, - ■ ■ ,Pa are connecting 
paths for the given quad. Figure [2] shows an example of this. It is easy to see 
that for any quad (Ci, C2, C3, C4) together with connecting paths pi, their union 
(Ui ^i) ^ (UiPi) contains as a minor. 




Figure 2. An example of a quad, with thickened connecting 
paths. 



We say a quad (Ci, C2, C3, C4) in a spatial graph F is double-linked (in Z or 
Z2) if lk(Ci, C3) ^ and lk(C2, C4) 7^ (in Z or Z2, respectively). 

It is easy to see that a spatial graph that contains a double-linked quad contains 
a double- linked Z?4 as a minor. Lemma [2] which we will prove later, says that the 
converse of this is also true, i.e.. 

If a spatial graph F contains a double- linked -D4 (in Z or Z2) as a 
minor, then it contains a double-linked quad (in Z or Z2, respec- 
tively) . 

Hence an embedding F of a graph G is a _D4-less embedding iff for every quad 
(Ci,C2,C3,C4) in F, lk(Ci,C3) = or lk(C2,C4) = 0. We translate this into 
systems of linear equations as follows. Let Fq be a fixed arbitrary embedding of 
G. Then any embedding F of G can be obtained from Fq by a sequence of crossing 
changes together with ambient isotopy in . Let ei, • • • , e„ denote the edges of 
G. For each pair of nonadjacent edges and e^ , let Xij denote the total number of 
crossing changes between and Cj in the given sequence of crossing changes that 
takes Fq to F. (Crossing changes between adjacent edges are ignored since they do 
not affect any linking numbers.) Note that Xij is the algebraic sum of the crossing 
changes between and e^-, where we give each edge an arbitrary fixed orientation, 
and associate a -1-1 to each crossing change that changes a left-handed crossing into 
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a right-handed one, and a —1 for the opposite case (see Figure [s]). (Alternatively, 
one can think of each crossing change as adding a full positive or negative twist 
between the two edges.) 




Figure 3. Changing a left-handed crossing (a) into a right-handed 
crossing (b) is denoted by -1-1; the reverse is denoted by —1. 

For each pair of disjoint cycles Ci and Cj in G, denote the linking number of Ci 
with Cj in Fq and F by lk(Ci,Cj, Fq) and lk(Ci,Cj,F), respectively. Then 

lk(C, , Cj , F) = lk(Q ,C,,To)+ 

It follows that C has a Z34-lcss embedding F iff for each quad (Ci, C2, C3, C4), at 
least one of the following equations is satisfied: 

(1) lk(Ci,C3,Fo)+ "^ki^O 
or 

(2) lk(C2,C4,Fo)+ Y "^ki^O 

where the linking numbers lk(Ci, Ci+2, Fq) are calculated numerically by the algo- 
rithm. Thus, we have: 

Proposition. Suppose a graph C has Q quads. Choosing one of equations ^ and 
for each quad gives systems of linear equations, with each system consisting 
of Q equations, such that G has a D^-less embedding (in TL or TLi) iff at least one 
of these systems of equations has a solution (in Z or Z2, resp.). 

Hence, if G has no Z)4-less embedding, the algorithm would a priori need to check 
every one of these 2*3 systems of equations to reach that conclusion. For the graphs 
we worked with, the number Q of quads typically ranged from a few thousand to 
hundreds of thousands, making 2^ much too large to allow this algorithm to be of 
practical use. Before we describe the shortcuts we used to speed up the algorithm, 
we discuss two more issues: (i) integer solutions versus Z2 solutions, and (ii) what 
the solutions to the equations tell us. 

Solving a system of equations in Z, as opposed to in Z2 or in M, requires different 
methods and can be much slower. So our algorithm first tries to solve each system 
of equations in Z2. If there is no solution in Z2, we know that there is no solution 
in Z either; so the algorithm moves on to the "next" system of equations. However, 
if a system of equations turns out to have a solution in Z2, then the algorithm 
determines whether the system also has a solution in K. If it does, and this solution 
happens to be in Z, the algorithm outputs this integer solution and reports that 
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G has a _D4-less embedding in Z. If the system has no real solution or if a real 
solution is found that is not in Z, the algorithm moves on to the "next" system 
of equations. In the end, if no Z2 solutions were found for any of the systems, 
the algorithm reports that the graph is IK. If at least one Z2 solution was found 
but never any integer solutions (we do not recall coming across such a case in our 
trials), we only conclude that G has a _D4-less embedding in Z2, (although it may 
also have one in Z). 

Now let us discuss the case when the algorithm finds an integer solution for one 
of the systems of equations. This solution tells us how many crossing changes to 
make between each pair of nonadjacent edges in order to go from Tq to a I?4-less 
embedding F. These crossing changes can be implemented in different orders and 
different ways up to isotopy, thus yielding different embeddings of G (Figure [4]). 
And all of these embeddings will be I?4-less since linking numbers depend only on 
the total numbers of crossings between nonadjacent pairs of edges. Of course, some 
or all of these Z34-less embeddings might contain nontrivial knots. However, every 
time we tried to, we succeeded in manually finding a knotless embedding by using 
these prescribed crossing changes as a "rough guide" . Our manual process was by 
no means algorithmic; it involved making choices in how to implement each crossing 
change. Making successful choices depended on both experience and luck! 





Figure 4. Two different ways (among infinitely many) that a "+1 
crossing change" (i.e., a full positive twist) can be introduced be- 
tween the oriented edges ab and cd. 
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3. The Shortcuts 

All of the work in this section applies to both Z and TL-i . With this understood, 
we will stop repeatedly saying "in Z or Z2". 

Let G be a graph. For each quad q = {C\,C2,Cz,Ci) in G, let (/(1, 3) and 
q(2, 4) denote the equations one gets in terms of the variables Xij by setting each 
of lk(Gi,G3) and lk(G2,G4), respectively, equal to zero. Let EQ be the set of all 
equations one obtains in this way from all quads in G, i.e., 

EQ = {q{i, i + 2) : i = 1, 2, g' is a quad in G} 

Note that each equation A G EQ may come from more than one quad. Let's say 
two equations A,B& EQ are related if there is a quad q such that {A, B} = 
{(7(l,3),g(2,4)}. Then, in any £)4-less embedding of G, at least one equation in 
each pair of related equations must be satisfied. 

One of our shortcuts is as follows. For each equation A G EQ, we find the set 
R{A) of all equations that are related to A. (Note that by definition A R{A).) 
In each _D4-less embedding of G, either A is satisfied or every equation in R{A) is 
satisfied. So if R{A), considered as a system of equations, has no solution, then 
A must be satisfied in every D4-less embedding of G; in this case we label A as 
an indispensable equation. Wc check each equation in EQ for whether it is 
indispensable, and each time we find a new one, we add it to the collection of all 
indispensable equations found so far and then check whether this new collection, 
considered as a system of equations, has a solution. If the collection has no solution, 
we conclude that G has no Di-\ess embedding, and the algorithm stops. Even 
when the collection does have a solution, this shortcut can still greatly speed up 
the algorithm in either finding a _D4-less embedding or determining that there is 
none: from among the 2^ systems of equations, we work only on those systems 
that contain all the indispensable equations. 

Remark: Among the graphs that we worked with, wc came across graphs G of 
both of the following types: (1) G has no £'4-less embedding but its indispensable 
equations are satisfiable as a system; (2) G has no £>4-less embedding and its 
indispensable equations are not satisfiable as a system. 

The second shortcut that we used in the algorithm is as follows. Let EQ = 
{A\,--- ,Az] (where the number of equations, z, is at most twice the number 
quads). Wc represent each of the 2^ systems of equations by a string S of Os and 
Is of length z, where the zth digit of 5 is 1 or according to whether Ai is or is 
not in the given system of equations. 

Of course wc arc not interested in all strings of Os and Is of length z\ we say such 
a string is valid if for each pair of related equations A, , Aj , the ith or jth digit of 
5 is 1. The algorithm progresses through the set of all valid strings of Os and Is 
of length z from smallest to largest in lexicographical order. The shortcut consists 
of trying to skip over collections of consecutive strings each of which is either not 
valid or corresponds to a system of equations that can be predicted, as described 
below, to have no solution. 

Suppose we test and determine a system of equations to have no solution. Often, 
the rows of the coefficient matrix of such a system turn out to be linearly dependent. 
Wc find a linearly independent subset of the rows that span the rowspacc. This 
corresponds to a subsystem of equations that still has no solution. Thus, any string, 
valid or not, that contains Is corresponding to this subsystem should be skipped. 
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This allows the algorithm to speed up greatly by skipping over large numbers of 
strings. 

We now describe the third shortcut we use in the algorithm. Let C be a cycle in 
a graph G, and suppose there exist vertices v,w in C such that the edge vw is in G 
but not in G. Then we say the edge vw is a chord for G . If G has no chords, we say 
it is chordless. A quad is said to be chordless if all of its four cycles are chordless. 
The shortcut is to only find and work with chordless quads instead of all quads. 
Obviously, if a graph G has a double-linked chordless quad in every embedding, 
then G is IK. However, a spatial graph that contains a double-linked quad may 
a priori contain no double-linked chordless quad. So an embedding that has no 
double-linked chordless quads may not be a D4-less embedding. It appears though 
that such cases are rare; in every instance that our algorithm found an embedded 
graph with no double-linked chordless quads, whenever we tried to we succeeded in 
manually finding a knotless (and, hence, Z?4-less) embedding for that graph. 



4. Proofs of the Lemmas 

Suppose _ff is a minor of a graph G, with a given sequence S of edge contractions 
that turn a subgraph G' of G into H. We say a path p' inG' C G corresponds to a 
path p in H if the given sequence S of edge contractions turns p' into p. (Note that 
a path may be a cycle or just one vertex.) There may be different subgraphs G' of 
G with different sequences S of edge contractions through which one can obtain H 
as a minor of G; and correspondence between two paths certainly depends on the 
choice of G' and S. So whenever we say p' corresponds to p, it is to be understood 
that the correspondence is relative to the given G' and S. 

When working with simple graphs, i.e., graphs with no loops or double edges, 
it is a common convention that whenever an edge contraction results in a double 
edge, the two edges in the resulting double edge are "merged together" into one 
edge. However, in this paper we work with D4, which is a multigraph; hence we 
do not adopt this convention. In other words, we do allow loops, and we do not 
merge together double edges. An advantage of this is that if H is obtained from 
a subgraph of a graph G by a finite sequence of edge contractions, then clearly to 
each edge in H there corresponds a unique edge in G. As we will see in Lemma [l] 
below (which is used to prove Lemma [2|, the same is true for cycles (though not 
necessarily for all paths). 

Lemma 1. Suppose H is obtained from a subgraph of a graph G by a given sequence 
of edge contractions, where double edges are not merged together and loops are not 
deleted. Then to every cycle C in H there corresponds a unique cycle C' in G. 

Proof. We use induction on the number of contracted edges. Suppose H is obtained 
from a subgraph G' of G by contracting one edge cq in G'. For each edge e in 
let e' denote the unique edge in G' that corresponds to e. Let G be a cycle in 
H. Let E{G) denote the set of edges in G. Then exactly one of Uee£;(C) ^' 
Ue6S(c) ^' ^ {^0} is a cycle in G' C G; and this cycle becomes G after we contract 
Cq. The lemma now follows by induction. □ 

Lemma 2. // a spatial graph F contains a double-linked D4 (in Z or 1^2) as a 
minor, then it contains a double-linked quad (in or 'Li, respectively). 
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Proof. We use induction on the number of contracted edges. To prove the "base 
case" , we simply note that a double-linked clearly contains a double-linked quad. 

As our induction hypothesis, assume is a spatial graph that contains a double- 
linked quad and is obtained from F by contracting one edge eo C F into a vertex 
vq €z fl. We will show F contains a double-linked quad. 

Let (Ci, C2, C3, C4) be a double-linked quad in D, with connecting paths pi, • • • ,P4. 
By Lemma[Tj there exist unique cycles C[, - ■ ■ , C4 in F that correspond to Ci, • • • ,6*4. 
It is easy to check that lk(Cj', C'^^i) = lk(Ci, C^+i) (in Z or Z2). We now construct 
connecting paths • • • , <Z4 for {C[, • • • C4). (We use qi instead of to denote these 
paths because some of them may not correspond to the pi.) Fix i S {1, 2, 3, 4}. We 
have two cases. 

Case 1. Ci n Ci+i is nonempty and connected. Then one verifies easily that 
C- n Cj'_|_j too is nonempty and connected. So we let qi be any vertex in C- n C^'^^. 
Case 2. d n Ci+i — 0. We have two subcases: 

(a) vq is not in pi, or vq is in the interior of pi. Then, by an argument similar 
to the proof of Lemma [l] there is a unique path in F that corresponds to 
Pi. We let qi be this path. 

(b) Vq is in Pi and vq is not in the interior of pi. Hence vq is an endpoint of 
Pi and must lie in exactly one of Ci or C^+i. Without loss of generality, 
assume vq G Ci. Let p'i be the path whose edges correspond to the edges of 
Pi. Upi has an endpoint in C.^', we let qi — p^; otherwise we let qi — pJU{eo}. 

Now it is not difficult to see that (C( , Cj, C3, C4) and gi,--- ,(74 satisfy the 
definition of a double-linked quad in F. □ 

Acknowledgment. The authors thank Thomas Mattman of CSU Chico, and the 
referees, for helpful comments. 
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